o ■ 
o ■ 

^ Photon density of states for deformed surfaces 

o : 

^ ■ T Emigi 

' ^ Institut fiir Theoretische Physik, Universitat zii Koln, Ziilpicher Strafie 77, 

D-50937 Koln, Germany 

• E-mail: te@thp.uiii-koeln.de 

o ■ 

(D ; 

^ , Abstract. A new approach to the Hclmholtz spectrum for arbitrarily shaped 

^ K ' boundaries and a rather general class of boundary conditions is introduced. We derive 

■ the boundary induced change of the density of states in terms of the free Green's 

C/3 I function from which we obtain both perturbative and non-perturbative results for the 

Casimir interaction between deformed surfaces. As an example, we compute the lateral 
electrodynamic Casimir force between two corrugated surfaces over a wide parameter 
range. Universal behavior, fixed only by the largest wavelength component of the 
' O ' surface shape, is identified at large surface separations. This complements known 

^ , short distance expansions which are also reproduced. 
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^ ■ 1. Introduction 

^ '• Casimir interactions can be viewed as a consequence of a change in the photon spectrum 
Y '. under a relative displacement of the interacting objects. For two parallel ideal metallic 
! plates the change in the photon density of states (DOS) with the plate distance H is well- 
O '. known It is peaked close to the characteristic frequency uj c/H which provides 
\ the dominant contribution to the interaction. However, even for simple geometries 
! like a plate and a sphere no exact expression for the change in the DOS is available. 
^ \ The reason for that is the difficulty to compute the distribution of eigenfrequencies of 
5^ \ the Helmholtz wave equation in arbitrary geometries. This is inherently linked to the 
problem of chaotic (quantum) billiards 12]. For Casimir interactions the natural question 
arises to what extent the force characterizes the shape of the interacting objects and 
vice versa. There is a long debate about the existence of a repulsive Casimir interaction 
between ideal metals which was mainly initiated by the positivity of the Casimir energy 
of a single infinitesimal thin sphere of ideal metal So far, no example of a repulsive 
force between two separated ideal metallic bodies has been found. 

During the last decades a number of approximative methods have been developed 
to determine the Casimir interaction in more complicated geometries. These include 
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the proximity and pairwise additivity approximations semiclassical approaches [H], 
multiple scattering expansions |I] and, more recently, a simple to implement optical 
approach [Hj which is however limited to short distances. Despite the usefulness of these 
approaches in certain limits, reliable results for the entire range of energies (or distances) 
are not known in general. Here I shall present a formula for the density of states 
in terms of the free space Green's function which is evaluated at the (ideal metallic) 
surfaces only. This formula will be used to compute the DOS for two deformed surfaces 
both perturbatively in the deformations and also exactly with the aid of a previously 
developed numerical method for periodic surfaces |7j. As an application, we study 
the lateral force between periodic surfaces. It is found that the surface's shape can 
be deduced from the force at short distances whereas at large separations universality 
prevails. 



2. Trace formula for the photon spectrum 

For geometries with a translational invariant direction the electromagnetic field splits 
into transverse electric (TE) and magnetic (TM) modes. Both modes are described by a 
scalar field which fulfills either Dirichlet or Neumann boundary conditions on the ideal 
metal surfaces. Hence, we consider in the following the scalar field Helmholtz equation 

(V2 + A;2)0(x) = O. (1) 

The relevant quantity for the Casimir interaction is the change 6p{k) = p{k) — poo{k) 
in the DOS caused by moving the surfaces from infinity to the finite distance for which 
one is interested in the Casimir energy. Hence Poo(^) is the DOS for infinitely separated 
surfaces. The relevant finite part of the zero-point energy is then given by the sum of 
kSp{k) over k. It is more convenient to work in Euclidean space by performing a Wick 
rotation to imaginary frequencies. With Sp{qo) = —6p{iqo) the Casimir energy becomes 



he f°° 

^ = Y J ^(lo(lo^p{(lo) ■ (2) 
In order to derive the DOS we start with the 4-dimensional Euclidean action 

S=\j d'X{V<pf. (3) 

Using a path-integral formulation with delta-functions enforcing the boundary 
conditions on the surfaces Sa 0, one obtains the two-point correlation function in 
the presence of boundaries p. With the free space Green's function Go(x, x';go) = 
e~'°l''~'''l/(47r|x — x'l), the modified correlations G(x, x';go) = (0go(^)</'-go(^')) the 
presence of boundaries is given by 

G'(x,x';go) -G'o(x,x;go) = (iuc/u'G'o(x, s„(u);go) 

a/3 

X M^liu, u';go) G'o(x', s^(u');go) , (4) 
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where 7W~^(u, u';go) is the functional inverse of A1q/3(u, u'; go) = G'o(Sa(u), s^(u'); go) 
with the 3D vectors s„(u) denoting the positions of surface Sa as a function of the 2D 
coordinate u. This result applies to Dirichlet boundary conditions. For the Neumann 
case Gq has to be replaced by 9n^(u)Go(x, SQ(u);go) on the right-hand side of Eq. (jH) 
and Maf3{u,u';qo) = dn^(^)dnf,{u')Go{sa{u),Sf3{u'); qo) with dn„ the normal derivative 
at the surface 5*^. Next we make use of the property that the trace of Green's function 
yields the DOS. In the Euclidean formulation this relation reads 



Here the spatial integration extends over all inner and other regions of the surfaces. Since 
Eq. (j3]) holds in any spatial region, we obtain the change in the DOS by integrating the 
right-hand side of that equation with x = x' running over the entire 3D space. The 
integration is easy to perform in momentum space since x and x' occur as arguments 
of the free Green's function in Eq. Using for the Fourier transformed propagator 
9/dqo{qQ + g^)~^ = — 2go(go + g^)^^, one obtains after x-integration both for 7V1 and 
AioD the final trace formula 



The trace runs over the 2D coordinates u and the surface indices a. An advantage of 
this formula is its independence of surface self-energies which diverge in the limit of ideal 
metals but are irrelevant for the interaction. It is important to note that this formula 
can be applied in any spatial dimension and also to other types of boundary conditions 
which can be implemented via delta-functions in the path integral jlOj . 

3. Density of states 

In this section we will apply the trace formula of Eq. (jHI) to two deformed but on 
average parallel plates. In the first part we study (in amplitude and curvature) small 
deformations of general shape (corresponding to the large distance limit) by perturbation 
theory. Non-perturbative numerical results for the DOS of periodically corrugated 
surfaces with edges at closer distances are presented in the second part. 

3.1. Perturbation theory for arbitrarily deformed surfaces 

Let us consider two parallel surfaces with small deformations which are parametrized 
by Si(u) = (u, /ii(u)) and S2(u) = {u, H + h2{u)) where the deformation fields ha are 
small compared to the mean distance H. We are interested in the deformation induced 
modification of the flat surface DOS 




(5) 



l_d_ 

Txdqo 



i:i\n{MM^) . 



(6) 



^Po(go) = f log(l-e-2''"^)/l 



(7) 



Photon density of states for deformed surfaces 



4 



where A is the surface area. Expanding the kernel Ai = Aio + SAi with A4q the kernel 
for flat surfaces and 6A4 the correction due to the deformations, one finds for the change 
of the DOS 

SpiQo) = --Tr [logA^o + log(l + Mo^6M)] - ''H oo" , (8) 
vr 

where the contribution for H oo has to be subtracted. Expanding the logarithm to 
second order in the deformations /^^(u), the change of the DOS can be written as 

,3 



+ / Uy2Kii{qo;u-W)[h^in)-h^iW)]' 



+ ^ i^x (go; u - u)ha{u)h(j{u) > 



(9) 

where Spo is given by Eq. (|7j). The kernels of the non-local contributions can be expressed 
in terms of the free Green's function. In the following we will focus due to limitations in 
space on Dirichlet conditions; the Neumann case has the form of Eq. Q with however 
different kernels. Defining the series 



^evenWo, " 2^ J^;^ 



n=l 
oo 



Rodd{qo;u) = ^ (2n-lY ^ ^ 



n=l 

the kernels can be written as 



^ll(?o;u) = — (9go 

TT 



Kx{qo;u) = - -(9go [dHRodd{qo;u)Y ■ (11) 



u 

\o-HRodd{qo;uy^ 

TC 

The contributions to the kernels can be viewed as multiple scattering paths between 
flat surface positions. The kernel K\\ has only contributions from an even number of 
paths crossing the gap, connecting two positions on the same surface, see Fig. [H^a). The 
first factor ~ OuGq in K\\ is represented as a non-crossing paths. An odd number of 
crossing paths yields the kernel K^, see Fig. Hl^b). Each path is associated with the free 
Green's function which is reduced by a factor which is the inverse squared number of 
gap crossings, cf. Eq. (fTU)). 

3.2. Non-perturbative treatment of periodic surfaces 

The above perturbative treatment is limited to small and smooth deformations. 
However, the trace formula turns out to be useful also in the case of surfaces with 
edges, and at large deformation amplitudes. Recently, a numerical tool has been 
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developed to evaluate the Casimir force for periodic surface deformations [7]. Here 
we shall demonstrate that same approach can be also applied to the DOS. Specifically, 
we consider uni-axial rectangular surface profiles shown in Fig. |2fa). In this section we 
restrict the analysis to no lateral shift, 6 = 0. By Fourier transforming 7Wq,/}(u, u'; go) 
with respect to u and u' one can employ the periodicity of the surface profiles (along 
the Xi-direction) which leads to the representation 

(p, q; go) = (27r)25 (p2 + ^2) 

(12) 



00 

X 

m=—oo 



J2 5iPi + qi + 27rm/X) N^^iq,,q2;qo) 



which defines the 2x2 matrices Nm that can be computed analytically and are given 
in [7j for the present geometry. j\4 can be transformed to block-diagonal form so that 
each block couples only waves whose momenta differ from a given Bloch momentum by 
integer multiples of l-njX. This decomposition can be used to define the function 

^(gi, qi\ go) = tr(i?-ia,„i? - B-^d,,B^) (13) 

where the trace runs over all discrete indices of the matrix B'^l^{qi, q^, go) = A^fcf/(gi + 
27r//A, ^2; go), and is the analog of A^oo, i-e., B for H ^ 00. For practical 
computations we take k, I = —M, . . . ,M with an integer cutoff M which should be 
taken to infinity. The change of the DOS can then be expressed as 

^pilo) = -^j dq2 j dqig{qi,q2;qo) . (14) 

The following results for electromagnetic fluctuations are based on the application of this 
formula to Dirichlet (TM modes) and Neumann (TE modes) boundary conditions where 
the inversion of B, the trace in Eq. ()13j) and the integration are performed numerically. 

Fig. 12 displays the change of the DOS separately for TM and TE modes at the 
fixed distance H = 10a and A ranging from a to 100a. Plotted is also the analytic 
expression for flat plates at the same mean distance and at the reduced distance H — 2a. 
An expected feature of these results is that the dominant contribution to the change 
of the DOS comes from frequencies go ~ H^^. For small modulation lengths A -C if 
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Figure 2. Change of the density of states 5p{qo) per unit area for TM and TE modes 
with parameters H = 10a and A/a = 1 (a), 10 (b), 100 (c). The results were obtained 
with a cutoff of M = 13. (For scahng of the (ifqo)-axis, H = 10a.) 
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one can expect that for go ^ the modes cannot probe the narrow valleys of the 
surface, and the DOS should be well described by the result for fiat plates at a reduced 
distance H — 2a. While this is fully confirmed for TE modes, the magnetic modes show 
an unexplained deviation at sufficiently small qo, see Fig.lSfa). For A ~ if, Fig-Efb), no 
significant differences between the two types of modes are observed, and the change of 
the DOS is in between the result for flat surfaces at distances H and if — 2a, respectively. 
For large A ^ ii. Fig. Efc), there is a low density of edges, and one can expect the DOS 
for flat plates at distance if to be a good approximation. Such behavior is indeed 
observed with a small decrease for TM and a small increase for TE modes close to the 
peak. 

4. Universality of the lateral Casimir force 

Finally, we use the previous results to study the Casimir force for the geometry in 
Fig. Efa) which emerges from the variation of the zero-point energy with the shift b. 
The energy can be obtained from the DOS via Eq. Without giving further details 
we note that the lateral force Fi = —dS/db can be obtained by computing directly the 
derivative of the matrices N^^ in Eq. ()12|) with respect to b. The results are shown 
as a function of the surface distance for a flxed b = A/4 in Fig. IHfc) and at a flxed 
H = 10a as a function of b in Fig. Efd). Let us compare to approximative methods 
for small H. The proximity force approximation |3] yields the lateral force i^ppA = 
A[2So{H)-So{H-2a)-So{H + 2a)]/XfoTO < b < A/2 with ^:o(ii) = - {1:^ 720) he / H^, 
which changes sign at b = A/2 discontinuously. The pairwise summation (PWS) of 
Casimir- Polder potentials U{r) = — {71/ 24) he/ r"^ is often naively applied to metals, using 
an "corrected" amplitude as to reproduce the correct result for flat ideal metal plates 
For the present geometry i^pws can be computed by a simple numerical integration. The 
two approximations are also shown in Fig. ini^c),(d). For small 6, both approximations 
agree and match our results. Beyond 6 > A/20 the PFA starts to fail since it does not 
capture the exponential decay of Fiat for increasing 6. The PWS approach has a slightly 
larger validity range [cf. Fig|2fb)] and reproduces the exponential decay. However it 
deviates by at least one order of magnitude from i^iat for 6 > 2.5A. 

At large 5 > A the lateral force shows universal behavior in the form of approaching 
a limiting form which is independent of the detailed shape of the surface corrugation. 
The universal limit is given by the lateral force for a sinusoidally shaped surface (with 
amplitude Oq and wavelength A) (TT], 



which is expected to hold for oq <^ A. Using the lowest harmonic of the present 
rectangular geometry, corresponding to ao = Aa/n, we flnd indeed excellent agreement 
between Fpt and our results for the geometry of Fig|2fa) for large S > X, see FiglHfc). 

For the dependence on b, see Fig. Ol^d), three regimes can be identifled. For A ^ if , 
the force proflle resembles almost the rectangular shape of the surfaces, and the PWS 



pt - 




(15) 
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Figure 3. (a) Two plate geometry, (b) Approximate validity ranges of PFA and PWS, 
and the universal parameter range. The vertical bar marks the parameter range for 
A/a in (d). (c) Lateral force (in units of the normal force Fo{d) between fiat surfaces). 
Shown are also the PFA (dash-dotted), PWS (dashed) and the universal limit Fpt 
(dotted), (d) Shape dependence of the lateral force. The dashed and the dotted curves 
are the PWS and the full perturbative result of JI], respectively. 



approach yields consistent results. For decreasing A, yet larger than H, the force profile 
becomes asymmetric with respect to 6 = A/4 and more peaked, signaling the crossover 
to the universal regime for X < H where the force profile becomes sinusoidal. We note 
that the PWS approach fails to predict the asymmetry of the force profile, and the PFA 
even predicts no variation for < 6 < A/2. 
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